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SUMMARY 


In  active  sonar,  large  time-bandwidth  product  echoes  from  submerged  objects  are 
oftv  .1  identified  through  the  use  of  a  matched  filter,  or  replica  correlator.  The  theoretical 
char.'.cteristics  of  replica  correlation  (RC)  are  well  known,  but  are  predicated  upon  the  echo 
being  identical  to  the  transmitted  pulse,  except  for  the  presence  of  white  Gaussian  noise. 
When  the  ocean  and  target  combine  to  filter  the  incident  pulse  in  a  random  manner,  the 
replica  may  no  longer  be  a  match  for  the  echo. 

Echo-Echo  Correlation  (EEC)  is  a  method  for  e.'.ploiting  the  (assumed)  stationarity 
of  lire  medium  and  target  to  allow  each  echo  in  a  train  of  three  or  more  closely  spaced 
echoes  to  be  used  as  a  match  for  the  following  echo.  This  results  in  a  correlation  technique 
that  is  insensitive  to  “echo  splitting,”  target  dopplcr,  and  pulse  type.  The  sampled  output 
of  the  correlator  may  be  reduced  by  an  almost  arbitrary  amount,  thereby  reducing  the  false 
alarm  rate,  but  without  affecting  the  probability  of  detection. 
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I 


INTRODUCTION 


The  problem  of  detecting  the  presence  of  an  underwater  target  began  to  receive 
attention  during  World  War  II  as  the  danger  from  submarines  and  subsurface  mines  began  to 
increase.  In  subsequent  years  the  problem  has  expanded  to  include  avoidance  of  submerged 
navigational  obstacles  such  as  wrecks  and  sea  ice.  The  primary  means  now  in  use  to  detect 
these  submerged  dangers  is  sonar,  which  functions  either  by  “listening"  (passive  reception) 
to  the  emitted  sound  of  a  target  or  by  actively  transmitting  sound  toward  the  target,  then 
listening  for  the  “echo,”  which  is  the  reflection  of  the  transmitted  energy  from  the  target. 
The  echo  will  return  to  tiie  receiver  distorted  to  some  degree  by  unwanted  energy,  which  is 
called  noise.  The  receiver  system  must  contain  some  form  of  decision-making  process  which 
analyzes  the  returning  energy  with  the  goal  of  identifying,  or  detecting,  the  echo  in  the 
presence  of  the  noise. 

An  echo  returned  from  a  perfect  reflector  will  be  identical  to  the  transmitted  pulse 
p(t),  provided  the  medium  through  which  the  energy  travels  does  not  distort  the  wave.  In 
the  actual  sonar  environment,  the  target  seldom  acts  as  a  perfect  reflector,  and  the  medium 
nearly  always  modifies  the  wave  in  several  ways.  A  transmitted  narrowband  pulse  of  sound 
will  be  characterized  by:  its  pulse  length  t,  its  center  frequency  fi,  its  bandwidth  B 
(B/J2<^1 ),  its  initial  phase  i//,  and  by  its  power.  The  parameters  may  be  modified  by  the 
target  and  the  medium  in  such  a  way  that  statistical  estimation  of  each  becomes  necessary. 
When  the  range  and/or  the  relative  velocity  of  the  target  is  unknown,  the  time  of  arrival  of 
the  echo  may  also  have  to  be  estimated. 

When  the  echo  is  so  buried  in  noise  that  its  immediate  identification  is  difficult, 
signal  processing  is  required  to  extract  the  echo.  The  primary  form  of  signal  processing  now 
in  use  with  sonar,  for  narrowband  pulses  in  which  the  product  Br  is  much  greater  than  1 ,  is 
called  Replica  Correlation  (RC).  Wlien  the  spectral  density  of  the  noise  n(t),  is  also  a  nar¬ 
rowband  function  centered  on  fi,  the  received  time  function  will  be 

z(t)  =  p(t)  +  n(t) 

=  R41[Z(t)  exp(int)l  (1) 

=  Rill(P(t)e‘’^  +  N(t))exp(int)| 


where 


Z(t),  P(t),  and  N(t)  are  the  complex  envelopes  (Ref.  (1 ))  of  z(t),  p(t),  and  n(t), 
respectively. 


The  observation  period  of  v(t)  is  T>r.  Assume  the  data  are  sampled  according  to  the  Shan¬ 
non  sampling  rule,  at  h  samples  per  second.  Replica  correlation  is  defined  as 


hT 


r(X)  pa  +  A)z(j) 

j=0 

'  hT 


hT 


'M 


^  P*0 +X)Z(i) 

j=0 


(2) 


in  which  pO)t  a  replica  of  the  transmitted  pulse,  is  the  optimum  form  of  the  matched  filter 
for  white  noise.  The  time,  Xmax,  of  the  largest  amplitude  of  r(X)  is  taken  as  an  estimate  of 
the  time  of  arrival  of  the  echo. 

Replica  correlation  presupposes  that  the  echo  is  identical  to  the  transmitted  signal 
save  for  the  effects  of  noise.  The  medium  and  target  will  often  combine  to  filter  the  trans¬ 
mitted  signal  in  such  a  way  that  the  received  echo  no  longer  duplicates  the  replica.  In  gener¬ 
al,  it  is  not  possible  to  predict  the  characteristics  of  this  “external  filter,’’  and  means  are 
required  to  compensate  for  its  effects,  although  the  techniques  employed  with  replica  corre¬ 
lation  have  not  been  notably  successful  when  severe  distortion  of  the  pulse  occurs.  The 
remainder  of  this  paper  presents  a  technique  that  eliminates  the  problems  caused  by  the 
external  filter. 


ECHO-ECHO  CORRELATION 

If  the  external  filter  is  fairiy  stationary  over  a  period  of  several  pulse  lengths,  it  is 
possible  to  compensate  for  the  pulse  distortion  by  transmitting  three  or  more  closely  spaced 
pulses,  identifying  one  of  the  returning  echoes,  and  using  it  as  a  “replica”  for  the  remaining 
pulses.  Let  the  transmitted  signal  be  a  series  of  pulses 

M 

a(t)  =  y  p(t-mL)  (.3) 

m=l 


where  L  >  r  is  the  pulse  repetition  period, 
M  >  3  is  the  number  of  pulses, 
and 


ML<T. 


The  received  sampled  time  function  from  a  stationary  target  now  is  x(j),  where 
xU)  =  b0)  + n(j)  ,  0<j<hT 

where 


b(j) 


M 


2 

m=l 


q(j  -  m  L) 


(4) 


qO)  is  the  received  echo  from  one  transmitted  pulse,  p(j). 

At  any  time  jg,  begin  observing  x(j)  in  segments  each  hL  samples  l<^ng.  Let  xj^  be  tlie 
segment.  Define  Echo-Echo  Correlation  (liEC)  as 


(hL-lX|) 

>  x^G +  X)  Xj^+jG) 

\  /  . 

^0 

/  T  \ 

"hL-lXl 

1  -1X1 

XjO+XjXk+iO) 

. 

hL  <  X  <  hL 

where  X]^G)  is  the  complex  envelope  ofxjjO)- 

There  are  two  cases  to  be  considered.  First,  assume  that  each  individual  echo,  qgO)>  K  =  1 , 

2,  3,  is  wholly  contained  within  the  adjacent  segment,  Xk+^_]G).  Refer  to  Fig.  1,  where  the 
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Figure  1 

envelope  of  xG)  is  shown;  the  noise  is  suppressed.  Assume  that  the  first  echo,  qjC).  is  con¬ 
tained  in  segment  Xj^,  that  the  individual  echoes  qj^  are  identical,  although  they  may  no 
longer  resemble  pO).  and  that  the  noise  in  each  segment,  ni^(j),  is  statistically  independent 
i'rom  signal  and  from  noise  in  adjacent  segments.  Eq.  (5)  is  used  to  obtain  sample  correla¬ 
tions  of  adjacent  segments.  There  are  seven  input  segments,  which  provide  six  correlation 
segments.  Tliey  are 


E1c|^_2(X)I  ”  G"  ^  E  {n|^_2G  +  (j)}  =  0  (6) 

j=0 

where  E  is  the  expectation  over  the  noise  distribution, 

J  =  (hL-|\!)  , 

Elci^.l(X)l  =  j^  E{nk.i(j  +  Miqi(J)  +  nka)l}  =  0  (7) 

j=0 

.1 

E[C|^t.\)j  E{tqi(t  +  X)  +  nj^(t  +  X)!  Iq2(t)  +  nj,+  j(t)]  }dt 

j=0 

(8) 

=  J  ^  qi  tt  +  X)  q2(t)  dt 


iWJ,y 


J 

E[ck+i(X)l  =  J  ^  E{[q.^C X)  +  iik+l^  [1130  +  11^+20)]  } 

(9) 

J 

=  7  ^  92(i  +  X)q3(j) 

j=o 

E[Ck+2(X)]  “  ^  E {[q3(j  +  X)  +  11^+2^  ^ 

j=0 

L 

E[Ok+3(X)]  =7^  E  {11^430 +  X)  11^+40)}  =  0  .  (11) 

j=0 

Note  that  noise  correlated  with  noise  [Eqs.  (6),  (11)1  gives  zero  correlation,  us  does  eclio 
correlated  with  noise  [Eqs.  (7),  (10)] ,  while  echo  plus  noise  correlated  with  echo  plus  noise 
gives  the  autocon'clation  of  the  (;cho  [.Eqs.  (81,  (9)| .  Eacii  correlation  segment  is  2L  sec 
long. 


In  the  second,  and  more  general,  case,  each  echo  will  be  split  between  segments,  as 
shown  in  Fig.  2.  Note  that  x^O)  contains  a  portion  of  q],  while  x^+i  (j)  contains  portions 
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Figure  2 


of  both  q I  and  q2,  etc.  As  before  Eqs.  (6),  (7),  and  ( 1  1 )  sliow  E  lck_2l  ~ 

=  E[ck+3]  =  0.  However,  since  XkO)  contains  a  part  of  q  i(j)  which  is  identical  to  the 
similar  part  of  q2(j)  in  Xk+l(j).  Elckl  will  be  a  partial  autocorrelation.  Segment  Xk+i(j) 
contains  a  complete  echo  (although  the  parts  are  rearranged)  as  does  seginenl  Xk+ifi)-  Thus 
E[ck+i(X)]  follows  Eq,  (9),  and  produces  a  full  autocorrelation.  Elck+2(X)1  again  pro¬ 
duces  a  partial  autocorrelation.  Since  one  of  the  two  ca.ses  will  always  occur  for  a  fixed 
target,  Eq.  (5)  wilt  always  produce  a  full  autocorrelation  function  for  at  least  one  value  oi‘  k. 

With  both  RC  and  EEC,  the  passage  through  the  correlati)r  of  a  sufficiently  strong 
echo  will  result  in  a  “peak”  in  the  output  which  exceeds  in  amplitude  the  surrounding  back¬ 
ground  caused  by  noise  passing  througii  the  correlator.  When  no  piior  knowledge  of  the 
target’s  location  is  known,  the  output  of  the  replica  correlator  must  be  computed  for  every 
value  of  X  to  determine  which  r(X)  is  largest  in  magnitude.  For  EEC,  on  the  other  liand. 

Fig.  1  and  Eq.  (8)  show  that  the  peak  must  occur  at  the  center  of  the  output  segment,  X=0, 
for  a  fixed  target.  Thus  all  values  of  c^fX),  for  X=A0  and  all  k,  need  nol  be  computed.  The 
reduction  in  cominitation  ovei  RC  is  apparent. 


6 


When  the  target  moves  with  constant  velocity  relative  to  the  transmitter/receiver, 
the  effect  is  to  change  the  pulse  repetition  period,  L.  If  v  is  the  relative  target  velocity,  and 
V  the  speed  of  sound,  then  L'  is  the  echo  repetition  period. 


The  individual  echoes,  qg(j)>  “ire  each  frequency  shifted  relative  to  the  transmitted  pulse, 
but  not  relative  to  each  other.  Since  v  is  unknown,  it  must  be  estimated;  however,  since 
V  =  (L'/L  -  1 )  V/2,  L'  determines  v.  If  information  about  the  target’s  maximum  speed  is 
available,  then  the  difference  |L-L'|  <  5.  Thus  it  becomes  necessary  to  evaluate  cj^CX)  only 
for  |X1  <  6,  and  the  time  X^ax  ^  determine  (L-L'),  hence  L’  and  v. 


SIGNAL-TO-NOISE  RATIOS 

One  of  the  primary  measures  of  tlie  effectiveness  of  a  signal  processing  technique  is 
its  ability  to  improve  thesignai-to-noise  ratio  of  the  input  zO)'  Define  the  input  signal-to- 
noise  ratio  as 

Si  =  F/aJ  (13) 


where 

F  is  tlic  mean  square  amplitude  of  the  echo,  qO))  und  is  the  variance  of  n(j). 
The  output  signal-to-noisc  ratio  for  RC  [Eq.  (2)1  is 


Sq  =  {Elr(X,j,ax)l 


2BT(Si) 


(14) 


wliere 

2 

Oj.  is  the  variance  ol  r(X) 
and  for  EEC  [Eq.  (5)1  it  is 

So  =  (E[c|<(0)1j2/„2^^^ 

^  2BT(Si)2 


(15) 


■I 


I 

1 

1 

\ 


\  \ 


where 


2 

a  is  the  variance  of  Cl,  .-7  (X). 

%+2 

The  derivations  of  these  formulas  arc  given  in  Ref.  2. 

It  is  interesting  to  note  that,  because  the  EEC  peak  occurs  at  Xj^^gx  ~  0>  fo*"  ^  target 
whicli  is  fixed  relative  to  the  receiver,  the  Sq  predicted  by  Eq.  (15)  can  be  doubled.  Let 


(16) 


\(X)  =  ^lck(X)  +  Ck(-X)l  . 

Define  the  new  output  signal-to-noise  ratio  to  be 

^0  =  E(ck(0)] /a^  (17) 


where 


var(c|^4.2(X))  . 


Now 


E[ck(0)l  =  ^E[c|^(0)  +  ck(0)j  =  ElCk(O)]  , 

so  that  the  signal  part  of  Sq  is  unchanged  from  that  of  Sq.  The  noise  variance  is  reduced  by 
a  factor  of  two,  however,  since 

Oj  =  var  + 

=  [va.  C|^+2(X)  +  var  ci^+2(-X)) 

=  ^  [2  var  • 

Since,  for  X  0,  t;k+2(X)  and  <^g+2("^)  I'ave  independent  wliite  Claussian  noise,  one  has 

(19) 

In  the  remainder  of  the  paper,  the  use  of  cg(X)  will  be  referred  to  as  tlie  enhanced  HliC 
technique. 


FALSE  ALARMS 


A  processor  decides  that  an  echo  is  present  when  the  output  is  suffieiently  high. 

The  decision  level  D  is  often  determined  as  a  function  of  one  or  more  of  the  noise  statistics, 
so  that,  when  noise  alone  enters  a  correlator,  the  possibility  exists  that  the  output  will  ex¬ 
ceed  the  threshold.  When  this  occurs,  the  results  are  interpreted  as  a  false  alarm. 


For  Br^  1 ,  the  Centrid  Limit  Theorem  (Ref.  3)  ensures  tliat  the  jirohability  density 
function  (pdf)  of  the  correlation  of  the  noise  will  be  Gaussian,  both  for  RC  1  Eq.  (2)1  and 
EEC  [Eq.  (5),  (16)]  However,  in  the  estimation  of  sonar  parameters,  it  is  often  tlie  case 
that  the  phase,  i//,  and  carrier  frequency.  J2,  are  unimportant.  Since  the  time  functions  ol 
interest  are  narrowband,  the  phase  and  carrier  of  tlie  echo  envelope  can  be  eliminated  by 
computing  the  magnitude  of  the  complex  functions  [liq.  (2),  (5),  and  (16)],  Tims 


^(X)  =  lr(X)|  ,  forRC 
P)I^(X)  =  |c|^(X)|  ,  for  EEC 
aj^(X)  =  |C[.(X)|  for  enhanced  EEC. 


(2 


These  are  the  correlation  “envelopes,”  which  will  follow  the  Rayleigh  distribution  (Ref.  4) 
so  that  the  probability  of  a  false  alarm,  PFA,  is,  for  RC 


PFA(0>Dj) 


=  exp 


and  for  EEC’ 


PFA{pk>l)2) 


PFACai-  >  D3) 


wliere 


(2 


(2; 

(23 


Dj,  i=  1,2,3  is  the  respective  threshold 


1 

u“  ,  i  -  1 ,2  is  the  respective  noise  variance. 

i 

By  forming  a  rath)  of  Eqs.  (22)  and  (23),  and  setting  D3  =  D2,  one  has 


for 


,  D2>0  , 

which  shows  tiiat  using  a^fX)  in  place  of  P[j(t)  will  always  give  a  reduced  PF  A  for  a  given 
tlireshold  level.  On  the  other  hand,  if  PFA  in  Eqs.  (22)  and  (23)  are  equal,  then  D3 
=  D2/\/T,  so  tliat  the  same  result  is  obtained  with  a  lower  threshold  for  the  enhanced 
technique. 


Jt  was  previously  mentioned  that  one  need  compute  C[-(X)  and/or  c^CX),  hence  pi^(t) 
and/or  a^CA),  only  within  some  range,  -6  <  X  <  5.  By  not  computing  the  correlation  for 
|X|  >  6,  the  opportunity  for  noise  only  output  to  exceed  the  tlireshold  is  reduecd.  Let 


whicli  is  seen  to  be  a  ratio  between  the  necessary  output  duration,  26,  and  tlie  maximum 
duration,  2L.  The  probability  of  a  false  alarm  is  thus  reduced  by  the  factor  j. 


COMPUTER  SIMULATION 

One  method  for  checking  the  validity  of  assumptions  made  in  deriving  theoretical 
equations  is  to  make  measurements  on  computer  simulated  data  and  to  compare  these  re¬ 
sults  with  the  theoretical  predictions.  To  that  end,  a  puke  was  generated  and  added  to 
bandlimited  white  Gaussian  noise.  This  input  time  function  was  used  as  /(t)  from  Eq.  (.  1 ) 
and  was  pa.ssed  ihrough  botli  the  replica  correlatoi^and  the  Echo-i'cho  Correlator.  The 
output  was  observed,  and  the  signal-to-iioise  ratio  S^,  measin-ed.  When  Br>  1 ,  the  output  of 
the  correlator  will  be  Ciaussian,  so  that  Sq  may  be  computed  from 

So  =  (■*^)'  06) 

where 

F  is  the  peak  of  the  correlation  output  when  pulse  (plus  noise)  has  passed  ihrough 
the  correlator  Y  is  the  sample  mean  of  the  output  noise,  a  is  the  sample  varianee  of 
the  output  noise. 

Equation  (26)  compares  with  Eqs.  (14)  .md  ( 15)  when  Y  =  0.  When  the  envelopes  of  zero 
mean  Gaussian  correlator  outputs  (Eqs.  (20))  arc  copsidered,  a  more  appropriate  measure  to 
use  is  (Ref.  5) 


where 


(27) 


F  and  Y  are  as  before. 

When  Sj  is  computed  as  the  mean  square  pulse  amplitude  divided  by  the  noise  variance,  Eq. 
(27)  is  an  asymptotically  unbiased  estimate  of  as  predicted  by  Eq.  ( 14)  or  ( 15)  as 
appropriate. 


The  pdf  ol'  the  noise  output  can  be  estimated  by  computing  a  histogram  of  the  data 
and  comparing  it  with  a  theoretic:'  '  say  by  a  “chi-squared  Goodness  of  Fit  Test” 
(Ref.  6).  This  estimate  can  be  imp  ;  averaging  together  several  histograms  of  inde¬ 

pendent  data.  'The  mean,  Y,  and  vai  ,e,  o-,  of  the  data  are  estimated  as 
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MM 


Y  = 


MM 


I  V(X) 


X=1 


unci 


-2  1 

"  "  MM- 


MM 


X=1 


(Y(X)-Y) 


(28) 


(29) 


where  MM  >  ht  and  Y(X)  are  samples  of  0(X),  Pk(X),  or  Q;ii(X),  as  appropriate.  The  False 
Alarm  Rate  K  is  estimated  by  counting  the  number  of  times  in  the  observation  period,  T, 
that  the  noise  only  i)ortion  of  the  correlation  exceeds  a  threshold,  D.  This  threshold  is 
usually  a  multiple  of  the  output  sample  noise  standard  deviation,  a. 


The  waveforms  used  in  the  sea  trials  were  either  Linear  Period  Modulated  (LPM) 
pulses,  or  Pseudo  Random  Noise  (PRN)  pulses.  See  Reference  7  for  a  discussion  ot  LPM 
signals.  A  set  of  three  identically  generated  LPM  pulses  were  generated  in  the  computer  and 
imbedded  in  bandlimitcd  white  noise.  Tne  pulse  bandwidth  was  15-70  Hz,  and  the  pulse 
length,  T,  was  0.853  sec.  The  pulse  repetition  period,  L,  was  1 .333  sec.  The  noise  band¬ 
width,  to  the  3-dB  points,  was  approximately  15-70  Hz.  See  Figure  3,  which  represents  the 
case  in  which  each  echo  is  wholly  contained  in  a  segment,  .xk(t)  (Fig.  1).  Equation  (20)  was 
evaluated  from  these  data  with  results  as  shown  in  Figs.  4  and  5,  for  RC  and  EEC^ 
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Figure  3.  Computer  simulated  x(l)  IRef.  Kq.  (4)]  . 
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respectively.  Ten  trials  using  the  same  pulse  train  mixed  with  independent  noise  samples 
were  conducted  and  the  lesults  averaged.  The  input  signal-to-noise  ratio  Sj,  calculated  as  the 
mean  square  pulse  amplitude  divided  by  the  input  noise  variance,  was  2.85  dB,  and  the 
average  output  signal-to-noise  ratio  Sq,  calculated  by  Eq.  (27).  was  22.36  ±  0.90  dB  for 
Pl^fi),  and  25.14  ±  0.8  dB  for  The  predicted  Sq,  using  Eq.  (1 5),  was  23.1  dB.  A 
similar  test  was  conducted  for  a  case  of  arbitrarily  located  echoes  (with  respect  to  segment 
orientation  -  Fig.  2).  with  results  very  nearly  the  same  as  just  described.  Figure  6  shows  a 
plot  of  the  probability  of  false  alarms  exceeding  a  threshold,  D,  versus  the  threshold,  given 
as  fractions  of  the  output  mean,  -  Eq.  (28).  Note  that  the  measured  data  points  fit  the 
theoretical  curves  very  closely. 


0  .5  1.0  l.b  2.0  2,5  3.0  3.5  4.0  4.5  5. 


THRESHOLD  (FRACTION  OF  MEAN) 
Figure  6.  Probability  of  false  alarms. 


COMPARISON  OF  ECHO  DETECTABILITY 

The  relative  ability  of  two  signal  processing  techniques  to  extract  a  recognizable 
echo  from  noise  is  a  measure  of  the  efficiency  of  one  technique  versus  the  other.  The  fol¬ 
lowing  comparison  method  is  due  to  Brown  (Ref.  2,  pp  2-30),  with  modifications  to  suit 
the  more  restricted  nature  of  this  inquiry. 


Define  the  False  Alarm  Rate,  K,  to  be  the  number  of  false  alarms  which  occur  in 
some  unit  of  time,  and  assign  K  to  both  techniques.  Determine  the  Sj  necessary  for  each 
techniciue  to  attain  this  K  for  a  50%  probability  of  detection.  For  Bt  >  1 ,  an  acceptable 
estimate  of  K  is 


K  =  PFA  •  B  . 
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(30) 


Using  Eqs.  (24),  (25),  and  (26),  we  have 

K  =  PFA(0>Di)-Bi  =  PFA(p>D2)-B2  =  PFA(a>D3)*B2  (31) 


where 


Bj  is  i  =  1,2  bandwidth. 

Solving  Eq.  (32)  for  Dj,  we  obtain 

/  2  A  \l/2 

Di  =  (-2a;;|en(K)/B,) 

/  2  -  VI- 

^2  =  (-2a;^2^n(K)/B2) 

/  0  .  \l/2  ^ 

D3  =  I  -4o‘^  Sin  (K)/B2  1  =>/2D2  .  (32) 

The  mean,  Yj^  and  the  variance  uj^,  of  the  Rayleigh  distributed  noise  out  of  the  correlalor/ 
detector  is  related  to  the  variance,  0“,  of  the  input  noise,  n(l)  [Eq.  (1)1  by  the  I'ollowing 
(Ret'.  5):  " 

Yr  =  a„v/7r/"2,aR  =  ^  ^n  =  "n  • 

An  “effective  signal-lo-noise  ratio,”  dj,  is  given  for  each  teciuiique  as 
dj  =  2(j  log  [(Uj- XRi)/aRil 

=  20  log  [(Dj  -  a„ixA/2)/|3o„j)  .  (34) 


Equation  (34)  is  equivalent  to  Eq.  (26).  In  order  to  determine  Sj  from  Eq.  (14)  or 
(15),  it  is  first  necessary  to  compute  Sq  from  Eq.  (27),  wliich  is  not  linearly  related  to  d|. 
The  non-linearity  may  be  removed  by  application  of  (modified)  Brown’s  curves,  lit  us  trans¬ 
forming  dj  to  an  estimate  of  Sq  comparable  to  that  obtained  from  E(|.  (27). 


As  a  numerical  example,  let  0,^ |  =  o„2  =  1 .0,  B  [  =  B2  =  B  =  400  llz,  r  =  1  sec,  so  1  0 
log  (2Br)  =  29  dB,  and  let  K  =  10"2  false  alarms  per  sec.  Thus,  from  Eqs.  (22),  (23)  and 
(24) 

Dj  =  D2  =  4.61 
D3  =  6.51 

and  from  Eq.  (34) 

dj  =  ds  =  d3  =  14,21  dB  . 

From  Fig.  7  we  obtain,  for  RC,  Sj  =  12-29  =  -17  dB,  while  from  Fig.  8,  for  EEC,  Sj  =  6.5- 
14.5  =  -8.5  dB,  and  for  the  enhanced  EliC  [  Eq.  (16)1,8;  =  5.5-1 4.5  =  -9.5  dB. 
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The  preceding  example  shows  RC  to  require  7.5  dB  to  8.5  dB  less  S;  to  attain  the 
given  K,  which  implies  that  RC  is  superior  to  EEC  as  a  detection  tool.  However,  this  exam¬ 
ple  required  two  assumptions.  First,  the  echo  was  not  distorted  by  the  external  filter  so 
that,  for  RC,  the  processing  gain  of  Eq.  ( 1 5)  holds.  Experimental  results  have  shown,  how¬ 
ever,  a  loss  of  up  to  IS  dB  from  the  anticipated  Sq  -  a  loss  that  would  not  occur  with  EEC. 
Second,  the  false  alarm  rate  assigned  for  RC  may  be  arbitrarily  reduced  for  EEC  simply  by 
not  computing  values  of  p|^(X)  or  agtX),  as  appropriate,  for  jXl  >  6.  Figure  9  shows  a  plot 
of  K  vs.  dj  for  different  values  of  7  [  Eq.  (25)1 .  To  complete  the  previous  example,  if 
7  =  0.01 ,  then  for  K  =  1 0"^,  the  EEC  technique  requires  that  dp  =  10.5  dB.  From  Fig.  8, 
the  required  S|  is  now  -11.5  dB,  or  an  equivalent  gain  in  detection  capability  of  3  dB. 


Figure  y.  False  alarm  rale  vs,  deleelioii  iiules. 


EXPERIMENTS  AT  SEA 


EXPERirdENTAL  CONDITIONS 

Tests  were  conducted  during  June  of  1977  at  a  site  located  approximately  at  46,5" N 
latitude,  12.8°W  latitude.  The  water  depth  was  apiiroximately  4150  metres,  and  the  sea 
state  was  a  moderate  three.  The  transmit/reeeive  platform  was  IIMS  MATAPAN,  a  emivert- 
ed  WW  11  battle  destroyer  with  a  side-mounted  planar  array  srmar.  The  submarine  target 
was  at  a  depth  of  1 00  ni.  The  majority  of  the  trials  occurred  during  parallel  course  (zero 
relative  velocity)  operations  with  a  target  aspect  angle  of  approximately  100  deg,  Tlie  sonar 
beam  deprcssion/elevation  angle  (D/E)  ranged  from  45  to  26  deg  relative  to  the  mean  sea 
surface.  The  acoustic  path  was  a  single  bottom  bounce,  and  llie  limiting  interference  about 
the  echo  was  reverberation  from  the  sea  surface  above  the  target  (sec  Fig.  10).  All 


transmissions  were  corrected  for  tlie  motions  of  the  platform,  but  the  received  signal  was 
not  so  corrected  prior  to  recording  on  magnetic  tape. 

SIGNAL  PROCESSING 


The  narrowband  signals  were  sampled  at  1024  samples/sec  —  approximately  2.5 
times  the  Nyquist  folding  frequency.  For  RC,  a  replica  of  the  transmitted  pulse  was  also 
sampled.  The  correlations  (Eq.  (20))  were  computed  via  a  Fast  Correlation  technique  (Ref. 
7),  and  the  results  plotted.  For  EEC,  the  output  segments  are  for  IM  <  6  =  10"^  seconds  of 
correlation.  The  input  segment  duration  was  the  pulse  repetition  period,  L  =  0.5  sec  or  0.6 
sec,  depending  on  the  lest,  giving  7  =  2  X  10“^,  or  1.67  X  10“^,  respectively. 

EXPERIMENTAL  RE  SULTS 

Although  the  object  of  this  experiment  was  to  compare  echo  detectability  using  RC 
and  EEC,  only  a  limited  amount  of  data  were  obtained,  and  the  errors  and  uncertainties  in 
the  data  made  quantitat'vc  measurements  impossible.  However,  a  qualitative  comparison  is 
possible  from  the  plots  that  follow.  These  plots  are  presented  in  groups,  with  the  first  plot 
in  each  group  being  a  representative  received  time  function*,  z(t)  [Eq.  (1)] ,  for  that  group 
of  transmissions.  Following  this  RTH  are  the  individual  RC  and  EEC  plots  t'or  each  trans¬ 
mission  (“ping”)  in  the  group.  The  groups  are  characterized  by  D/E  angle  and  pulse  type,  as 
indicated  in  Table  1 .  Each  RTH  has  a  heavy  horizontal  line  marked  “Echo  Time.”  Because 
the  target  was  equipped  with  a  delayed  acoustic  transponder,  it  is  possible  to  determine 
rather  precisely  the  expected  echo  time  from  the  target.  Each  c  rrelation  plot  has  a  similar 
line  marked  “Period  of  interest,”  indicating  the  anticipated  echo  eonelation  time. 

Figure  1  lb  shows  EEC  output  as  anticipated  from  theory.  There  is  a  single  strong 
peak  with  two  smaller  side  peaks  equally  spaced  on  either  side  of  it.  These  three  peaks 
occur  precisely  at  X  =  0  for  their  respective  segments.  Note  that  the  corresponding  RC  plot. 
Fig.  1  la,  does  not  show  an  acceptably  high  enough  echo  corrclat'on  to  reasonably  declare 
an  echo  to  be  present.  Most  of  the  remaining  plots  in  this  group  show  more  than  the  antici¬ 
pated  two  or  three  peaks  for  EEC,  which,  it  is  assumed,  is  accounted  for  by  some  unknown 
acoustic  paths.  The  nature  of  these  paths  is  beyond  the  scope  of  this  report.  It  is  apparent. 


*Tlie  uiceived  time  function  here  is  referred  to  as  a  Reverberation  Time  History  (RTH).  Since  the  back¬ 
ground  interference  for  tlicse  tests  is  dominated  by  acoustic  returns  from  environmental  scattcrers,  hence 
"reverberation.” 
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Tabic  1 .  Experimental  parameters. 


Figure 

D/E 

BW 

IT 

Ping  Numbers 

11 

26 

400 

LFM 

43X-442 

12 

26 

600 

U’M 

443-447 

13 

35 

500 

I’RN 

224-225 

14 

35 

400 

U>M 

210-223 

IS 

40 

400 

LPM 

147-155 

16 

45 

! 

400 

LPM 

_ 1 

82-87 

Where;  D/li  in  degrees,  UW  is  pulse  bandwidth  in  Hz,  and  PT  is  pulse  type. 


however,  that  at  this  D/E,  EEC  provides  a  much  tnore  .satisl'actory  target  itidicatioti  than 
docs  RC,  which  often  has  the  echo  correlations  nearly  hitried  in  the  backgroiitid. 

At  steeper  D/E  angles  (shorter  ranges),  the  KC  plots  begiti  to  improve  see  h'ig.  14e. 
However,  the  EEC  plots  are  still  superior  in  Sq  see  Fig.  14d.  Several  of  the  EEC  plots, 
such  as  14f  and  14h,  were  suppressed  by  the  large  spike  at  the  left  hand  side,  but  the  Sq  is. 
in  fact,  higher  than  that  for  the  corresponding  RC, 

At  the  steeper  D/E  angles  of  35  and  45  deg.  the  high  output  from  EEC  in  the  first 
seven  to  fifteen  segments  indicates  correlated  echoes  from  the  bottom.  As  an  example,  see 
Fig.  1  5h.  These  returns  may  serve  to  suppress  the  scale  for  the  remainder  of  the  plot,  but 
the  geometry  of  Fig.  10  shows  them  to  occur  well  before  any  possible  target  echoes,  so  that 
they  do  not,  in  fact,  affect  the  signal-*o-noise  ratio  of  the  target  echo. 


Figure  16,  which  shows  the  results  for  the  steepest  D/E  of  45  deg,  not  only  has  more 
correlation  peaks  for  both  RC  and  EEC  than  anticiimted,  but  tlie  liEC  peaks  are  not  cen¬ 
tered  at  X  =  0  for  their  respective  segments.  While  a  definitive  reason  for  either  anomaly  is 
not  apparent,  the  off-center  peaks  were  most  likely  caused  by  the  vertical  acceleration  of 
the  receive  platform  during  the  reception  time  of  tlie  target  echoes.  F'or  the  nu)st  shallow 
D/E  of  26  deg,  any  vortical  acceleration  of  the  receive  platform  would  have  produced  only  a 
minor  projection  along  the  beam  axis,  and,  as  shown  in  Fig.  1  1 ,  the  ireaks  are  constrained  to 
X  =  0. 


The  enhancement  process  iliscussed  earlier  (Etp  (20))  was  not  used  on  these  data 
because  of  the  peaks  not  being  constrained  to  segment  center  (X  =  0). 

In  order  to  compare  empirical  results  with  theoretical  predictions,  one  reciuirement 
is  that  the  background  noise  statistics  be  the  same  as  those  used  in  the  predictions.  1  he 
various  RTH  plots  show  pronounced  variations  in  background  intensity,  so  that  (he  concept 
of  noise  stationarity  does  not  hold.  However,  it  is  pi)ssiblc  to  reduce  the  effects  of  a  slowly 
time  varying  background  by  the  number  of  methods.  One  might  pass  the  correlator/ 
detector  output  through  a  suitable  high  pass  filter,  or,  perhaps,  since  (he  correlation  may  he 
accomplished  with  the  cross  spectrum,  by  dividing  the  cross  spectrum  by  the  auto  irower 
spectra  of  both  components  in  the  correlation,  ll  appears,  however,  that  (lie  EEC  plots  have 
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Figure  Ilf.  HEC (ping 445). 
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Figure  llg.  RC  (ping  446). 
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Figure  lii.  KC  (ping  447). 
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I’igure  llj.  HFC  (ping  447). 
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Figure  14g.  RC  (ping  222). 
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Figure  16a.  RC  (ping  82). 
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Figure  16b.  MFC  (ping  82). 
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a  rulhcr  characteristic  appearance  (see  Fig.  1 1,  for  example)  not  available  with  RC,  and  such 
a  “normalization”  would  destroy  this  appearance.  Further  work  is  required  on  this  point. 

Although  the  data  are  sparse,  the  EHC  results  appear  to  improve  with  respect  to  RC 
as  the  range  increases.  Compare  Fig.  1 1,  the  longest  range,  with  Fig.  16,  the  shortest.  This 
is  perhaps  a  result  of  the  increased  multipath  interference  with  range  degrading  RC. 

All  of  the  previous  discussion  presuppo.sed  identical  input  signals  to  (white)  noise 
ratios  as  the  basis  for  SNR  comparison  between  EEC  and  RC.  In  the  case  of  practical  sonar 
comparisons,  it  must  be  pointed  out  that  a  transmission  of  three  pulses  suitable  for  EEC 
requires  substantially  more  signal  power  transmission  than  does  the  .single  pulse  necessary 
i'or  RC  for  a  comparable  input  echo  to  reverberation  level,  Sp  The  EEC  pulse  train  causes  a 
three-fold  increase  in  background  reverberation  for  each  individual  returning  echo.  Against 
this  loss  must  he  measured  the  gain  inherent  in  multipath  recombination  that  is  characteris¬ 
tic  of  EEC,  as  well  as  the  reduction  in  false  alarm  rate. 


CONCLUSIONS 


Echo-Echo  Correlation  has  a  number  of  characteristics  that  distinguish  it  from  repli¬ 
ca  correlation  and  which  suggest  that  EEC  may  be  a  superior  detection  tool  under  a  variety 
of  environmental  and  ojrerational  situations.  A  summary  of  the  characteristics  is  given 
below : 

LOCATION  OF  TARGET  RETUI^N 

The  peak  echo  correlation  will  always  be  found  at  a  known  or  estimable  location  in 
the  correlator  output:  at  .segment  center  for  the  i'ixed  target  case  or  within  a  short  time 
span  about  the  segment  center  for  the  moving  target  condition.  Such  a  |)riori  knowledge  is 
not  available  I'or  replica  correlation, 

PLATFORM  STABILITY 

Sonar  platform  motion  will  shift  the  output  echo  peak  about  the  segment  center  in  a 
possibly  random  manner.  Without  a  stabilized,  or  motion  compensated  platform,  the  analy¬ 
sis  of  target  doppler  would  be  difficult. 

MULTIPATH  AND  ECHO  SPLITTING  LOSS 

EEC  is  relatively  insensitive  to  multiple  returns  of  both  types  provided  the  target  is 
not  undergoing  high  accelerations.  Replica  correlation  may  be  severely  degraded, 

TARGET  DOPPLER  SENSITIVITY 

EliC’  may  suffer  a  minor  loss  in  correlator  amplitude  caused  by  echo  truncation,  but 
is  insensitive  to  velocity-induced  frequency  distortion. 
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PULSE  TYPE  SENSITIVITY 


For  a  motion  compensated  platform,  EEC  is  wholly  insensitive  to  uulse  type,  As  an 
example,  doppler-distorted  PRN  is  handled  with  the  same  filters  as  any  for  other  pulse  type, 
but  the  replica  correlator  must  have  a  large  bank  of  doppler-shil  ted  replicas  to  handle  just 
this  case  alone. 

FALSE  ALARMS 

Since  the  peak  return  occurs  about  the  segment  center,  the  probability  of  detection 
is  independent  (within  doppler  constraints)  of  the  number  of  samples  displayed  in  the  out¬ 
put.  Thus,  only  a  small  portion  of  the  output  for  each  .segment  need  be  displayed,  thereby 
greatly  reducing  the  false  alarm  rate.  Replica  correlation  rciiuires  that  all  output  samples  he 
displayed,  with  a  much  higher  false  alarm  rate. 

SNR  ENHANCEMENT 

Knowing  the  precise  location  of  the  correlator  peak  allows  the  use  of  Ei).  (16), 
namely,  c|^(X)  =  cit(X)  +  ci^(-X)  where  Cj,^(X)  is  the  standard  EEC  output  in  a  segment, 

%(X)  will  have  a  3  dl)  gain  in  signal-to-noise  ratio  over  ci^(X). 

CHARACTERISTIC  APPEARANCE 

Replica  Correlation  in  a  reverberation  background  provides  distinguishable  detec¬ 
tions  only  in  terms  of  amplitude  above  background.  EEC  returns,  in  addition,  may  offer  a 
characteristic  apiiearance  that  is  a  useful  aid  in  detection. 
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